Magnetic Monopoles in Spin Ice 
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Electrically charged particles, such as the electron, are ubiquitous. By contrast, no elementary particles 
with a net magnetic charge have ever been observed, despite intensive and prolonged searches^. We pursue 
an alternative strategy, namely that of realising them not as elementary but rather as emergent particles, i.e., 
as manifestations of the correlations present in a strongly interacting many-body system. The most prominent 
examples of emergent quasiparticles are the ones with fractional electric charge e/3 in quantum Hall physics^. 
Here we show that magnetic monopoles do emerge in a class of exotic magnets known collectively as spin 
ice^i^: the dipole moment of the underlying electronic degrees of freedom fractionalises into monopoles. This 
enables us to account for a mysterious phase transition observed experimentally in spin ice in a magnetic field^'^, 
which is a liquid-gas transition of the magnetic monopoles. These monopoles can also be detected by other 
means, e.g., in an experiment modelled after the celebrated Stanford magnetic monopole search^. 
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Spin-ice materials are characterised by the presence of 
magnetic moments fli residing on the sites of a pyrochlore 
lattice (depicted in Fig. [T]i- These moments are constrained 
to point along their respective local Ising axes (the di- 
amond lattice bonds in Fig. [U, and they can be modelled 
as Ising spins fli — jiSi, where Si = ±1 and /i = \fli\. 
For the spin ice compounds discussed here, Dy2Ti207 and 
Ho2Ti207, the magnitude n of the magnetic moments equals 
approximately 10 Bohr magnetons (/i = lO/is). The thermo- 




FIG. 1: The pyrochlore and diamond lattices. The magnetic mo- 
ments in spin ice reside on the sites of the pyrochlore lattice, which 
consists of corner- sharing tetrahedra. These are at the same time the 
midpoints of the bonds of the diamond lattice (black) formed by the 
centres of the tetrahedra. The ratio of the lattice constant of the dia- 
mond and pyrochlore lattices is ad /a = \/3/2. The Ising axes are 
the local [111] directions, which point along the respective diamond 
lattice bonds. 



dynamic properties of these compounds are known to be de- 
scribed with good accuracy by an energy term that accounts 
for the nearest neighbour exchange and the long ranged dipo- 
lar interactions'*'^'"': 
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The distance between spins is 7-^ , and a ~ 3.54 A is the py- 
rochlore nearest-neighbour distance. D = fiQfj? /{Ana^) — 
1.41K is the coupling constant of the dipolar interaction. 

Spin ice was identified as a very unusual magnet when it 
was noted that it does not order to the lowest temperatures, T, 
even though it appeared to have ferromagnetic interactions^. 
Indeed, spin ice was found to have a residual entropy at low 
T^, which is well-approximated by the famous Pauling en- 
tropy for (water) ice, S ^ Sp = (1/2) log(3/2) per spin. 
Pauling's entropy measures the huge ground state degeneracy 
arising from the so-called ice rules. In the context of spin ice, 
its observation implies a macroscopically degenerate ground 
state manifold obeying the "ice rule" that two spins point into 
each vertex of the diamond lattice, and two out. 

We contend that excitations above this ground state man- 
ifold, i.e., defects that locally violate the ice rule, are mag- 
netic monopoles with the necessary long distance properties. 
From the perspective of the seemingly local physics of the ice 
rule, the emergence of monopoles would at first sight seem 
rather surprising. To demystify it we will probe deeper into 
how the long range magnetic interactions contained in Eq. [T] 
give rise to the ice rule in the first place. We then incorporate 
insights from recent progress in understanding the entropic 
physics of spin ice, and the physics of fractionalisation in high 
dimension s of which our monopoles will prove to 
be a classical instance. 

Consider a modest deformation of Eq. [T] loosely inspired 
by Nagle's old work'* on the 'unit model' description of water 
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FIG. 2: Mapping from dipoles to dumbbells. The dumbbell picture 
(c,d) is obtained by replacing each spin in (a,b) by a pair of oppo- 
site magnetic charges placed on the adjacent sites of the diamond 
lattice. Left: two neighboring tetrahedra obeying the ice rule, with 
two spins pointing in and two out, giving zero net charge on each 
site. Right: Inverting the shared spin generates a pair of magnetic 
monopoles (diamond sites with net magnetic charge). This config- 
uration has a higher net magnetic moment and it is favoured by an 
applied magnetic field oriented upward in the figure (corresponding 
to a [111] direction). Large panel: A pair or separated monopoles 
(large red and blue spheres). A chain of inverted dipoles ('Dirac 
string') between them is highlighted, and the magnetic field lines are 
sketched. 



ice: we replace the interaction energy of the magnetic dipoles 
living on pyrochlore sites by that of dumbbells consisting of 
equal and opposite magnetic charges that live at the ends of 
the diamond bonds (see Fig. |2]i. The two ways of assigning 
charges on each diamond bond reproduce the two orientations 
of the original dipole. Demanding that the dipole moment 
of the spin be reproduced quantitatively fixes the value of the 
charge at zLfi/ad, where the diamond lattice bond length ad — 

The energy of a configuration of dipoles is computed as the 
pairwise interaction energy of magnetic charges, given by the 



magnetic Coulomb law 

VM = f (2) 

where Qa denotes the total magnetic charge at site a in the di- 
amond lattice, and rap is the distance between two sites. The 
finite "self-energy" uo/2 is required to reproduce the net near- 
est neighbour interaction coiTectly. This expression - derived 
in detail in the annexed supplementary notes - is equivalent 
to the dipolar energy Eq.[T]up to corrections which are small 
everywhere, and vanish with distance at least as fast as 1 /r^. 

Consider first the ground states of the system. The total 
energy is minimised if each diamond lattice site is net neutral, 
i.e., we must orient the dumbbells so that Qa = on each site. 
But this is just the above-mentioned ice rule, as illustrated in 
Fig.|2] Thus, one of the most remarkable features of spin ice 
pops right out of the dumbbell model: the measured low-T en- 
tropy agrees with the Pauling entropy (which follows from the 
short-distance ice rules), even though the dipolar interactions 
are long-ranged. 

We now turn to the excited states. Naively, the most ele- 
mentary excitation involves inverting a single dipole / dumb- 
bell to generate a local net dipole moment 2/i. However, this 
is misleading in a crucial sense. The inverted dumbbell in fact 
coiTesponds to two adjacent sites with net magnetic charge 

Qa = ±qm = ±2^i/ad, (3) 

a nearest neighbour monopole-antimonopole pair. As shown 
in Fig. 12] the monopoles can be separated from one another 
without further violations of local neutrality by flipping a 
chain of adjacent dumbbells. A pair of monopoles sepa- 
rated by a distance r experiences a Coulombic interaction, 
—fioq^/ (ATTr), mediated by monopolar magnetic fields, see 
Fig.E] 

This interaction is genuinely magnetic, hence the presence 
of the vacuum permeability /io, and not 1/eo- It takes only a 
finite energy to separate the monopoles to infinity (i.e., they 
are deconfined), and thus they are the true elementary excita- 
tions of the system: the local dipolar excitation fractionalizes ! 

The reader may now feel that we have slipped the baby in 
with the bathwater - that in taking the pictures from the dumb- 
bell representation seriously we are somehow introducing 
monopoles where there were none to begin with. In general, it 
is of course well known that a string of dipoles aiTanged head 
to tail realises a monopole-antimonopole pair at its endsi^. 
However, in order to obtain deconfined monopoles, it is es- 
sential that the cost of creating such a string of dipoles remain 
bounded as its length grows, i.e., that the relevant string ten- 
sion vanish. This is evidently not true in a vacuum (such as 
that of the universe) where the growth of the string can only 
come at the cost of creating additional dipoles. Magnetic ma- 
terials, which come equipped with vacua (ground states) filled 
with magnetic dipoles, are more promising. However, even 
here a dipole string is not always a natural excitation, and 
when it is - e.g., in an ordered feiTomagnet - a string of in- 
verted dipoles is accompanied by costly domain walls along 
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FIG. 3: Monopole interaction. Comparison of the magnetic 
Coulomb energy — /iogm/(47rr) (Eq.[2] solid line) with a direct nu- 
merical evaluation of the monopole interaction energy in dipolar spin 
ice (Eq. [T] open circles), for a given spin-ice configuration (Fig.|2j, 
as a function of monopole separation. 

its length (except, as usual, for one-dimensional systems'^), 
causing the incipient monopoles to remain confined. 

The magic of spin ice arises from its exotic ground states. 
The ice rule can be viewed as requiring that two dipole strings 
enter and exit each site of the diamond lattice. In a typical spin 
ice ground state, there is a "soup" of such strings: many dipole 
strings - of arbitrary size and shape - can be identified which 
connect a given pair of sites. Inverting the dipoles along any 
one such string creates a monopole-antimonopole pair on the 
sites at its ends. The associated energy cost does not diverge 
with the length of the string, unlike the case of an ordered 
ferromagnet, as no domain walls are created along the string, 
and the monopoles are thus deconfined. 

Notice that we did not make reference to the famous Dirac 
condition'^ that the fundamental electric charge e and any 
magnetic charge q must exhibit the relationship eq ~ nh/fiQ 
whence any monopoles in our universe must be quantised in 
units of = h/fiQC. This follows from the monopole be- 
ing attached to a Dirac string which has to be unobservable.^^ 
By contrast, the string soup characteristic of spin ice at low 
temperature makes the strings energetically unimportant, al- 
though they are observable and thus not quantised. 

Indeed, the monopoles in spin ice have a magnitude 

2fi 

q-m = — 
ad 

= — 7, — qo 



where \c is the Compton wavelength for an electron, and 
a is the fine-structure constant. Amusingly, the charge of a 
monopole in spin ice can even be tuned continuously by ap- 
plying pressure, as this changes the value of /i/a^. 

The monopoles are sources and sinks of the magnetic field 
H, as is appropriate to the condensed matter setting. More 
precisely, like in other instances of fractionalisation-", one can 



define a 'smeared' magnetic charge 

p„(R)^|dVe-l'-'-^l'/«' V-H, (5) 

where V • H is the divergence of the magnetic field. For a 
monopole at the origin, separated by L ^ ^ a from any 
other monopoles, this gives Pm{Q) = ±?m- The form of the 
magnetic induction B is also monopolar, but with the impor- 
tant difference that a compensating flux travels along the (un- 
quantised) "Dirac string" of flipped dipoles created along with 
the monopole (see Fig. |2]i. 

Our magnetic monopoles would in principle show up in one 
of the best-known monopole searches, the Stanford experi- 
ment to detect fundamental magnetic monopoles from cosmic 
radiation. This experiment is based on the fact that a long- 
lived current is induced in a superconducting ring when a 
monopole passes through it**. One can easily check that the 
presence of the Dirac string of flipped dipoles is immaterial to 
the establishment of a current. 

The above observations are the central qualitative re- 
sults of our work: ice-rule-violating defects are deconfined 
monopoles of H, they exhibit a genuine magnetic Coulomb 
interaction, Eq. |2] and they produce Faraday electromotive 
forces in the same way as elementary monopoles would. 

We reemphasize that the ice rule alone does not permit a 
consistent treatment of the excited states of the physical prob- 
lem: crucially, the energetic interaction between our defects is 
absent altogether Also, in previous discussions of the purely 
ice-rule problem and related short ranged problem s "-^^''^ it 
has been noted that the defects do acquire a purely entropic 
Coulombic (i.e., 1/r) interaction that has a strength that van- 
ishes proportionally to T at low temperatures. This interaction 
will be present in addition to the magnetic Coulomb interac- 
tion discussed in this paper, and is clearly much smaller as 
T 0. Also, it will not he accompanied by a magnetic field, 
it will not renormalise the monopole charge, and it will not 
be felt by a stationary magnetic test particle embedded in the 
lattice but not attached to a lattice site. 

The most satisfactory way to demonstrate the presence of a 
monopole would be to measure the force on magnetic test par- 
ticles, say by a Rutherford scattering experiment or by clever 
nanotechnological means. Unfortunately, for lack of elemen- 
tary magnetic monopoles, one would have to choose dipoles 
as test particles, which significantly weakens such signatures. 

An alternative strategy is to look for consequences of the 
presence of magnetic monopoles in the collective behaviour of 
spin ice. This is most elegantly achieved by applying a mag- 
netic field in the [111] crystallographic direction. Such field 
acts as a (staggered) chemical potential, see Fig.|2] favouring 
the creation of monopoles of a given sign on either sublattice 
of the diamond lattice. 

We thus have a tuneable lattice gas of magnetic monopoles 
on the diamond lattice. The basic structure of the phase 
diagram as a function of magnetic field and temperature 
can be inferred from a series of papers by M. E. Fisher 
and collaborators^' in the context of ionic lattice gases and 
Coulombic criticality. At high T, there is no phase transition 
but a continuous crossover between the high- and low-density 
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phases as the chemical potential is varied. At low T, a first- 
order phase transition separates the two regimes. This tran- 
sition terminates in a critical point at {hc,Tc), not unlike the 
liquid-gas transition of water This serves as a useful diagnos- 
tic as the liquid-gas transition is absent for a nearest-neighbour 
spin-ice model, in which defects interact only entropically. In 
that case, it is known that there cannot be a first-order transi- 
tion in the limit of low T^^. 

To confirm this scenario, we have demonstrated by Monte 
Carlo simulations that the actual phase diagram of dipolar spin 
ice model has precisely this structure. In order to rule out the 
appearance of the liquid-gas transition due to effects intro- 
duced by the approximations leading to Eq. |2l we simulated 
directly the original dipolar spin ice model, Eq. [T] The result- 
ing phase diagram is depicted in Fig.|4] The critical endpoint 
is located around (T^, he) = (0.57 ± 0.06K, 0.86 ± 0.03T). 
The error bars are mainly due to finite-size effects, as the in- 
tensive nature of the simulations of long-range dipolar inter- 
actions prohibits simulating very large systems. 

This scenario is indeed observed experimentally in spin-ice 
materials6 7, and our results provide a natural explanation. 
Spin ice in a [111] magnetic field is a very interesting prob- 
lem which has already attracted considerable attention. The 
low-density phase of monopoles is known as kagome ice, a 
quasi two-dimensional phase with algebraic correlations and 
a reduced residual entropy^-^'^-^ . The high-density phase is an 
ordered state with maximal polarisation along the field direc- 
tion. Experimental results on the liquid-gas transition and its 
endpoint are also displayed in Fig.|4]for comparison. Our nu- 
merical results are in good qualitative agreement with both ex- 
periment and the analytic calculations of Ref . 2 1 . Our value of 
the critical field agrees with Ref. 6 within a few percent, less 
than the uncertainty due to demagnetisation effects^. How- 
ever, the experimental value of Tc is about a third lower than 
the numerical one, most likely due to farther-neighbour (ex- 
change) interaction terms, which despite being small can shift 
the location of a transition temperature considerably"'. 

The presence of a liquid-gas transition was noted to be 
very remarkable as there are few, if any, other experimentally 
known instances in localised spin systems'". No mechanism 
was known to account for this phenomenon, and our theory of 
magnetic monopoles fills this gap. 

The existence of magnetic monopoles in a condensed mat- 
ter system is, given their intellectual pedigree, exciting in 
itself. (The reader familiar with the beautiful work on the 
anomalous Hall effect will notice that the monopoles appear- 
ing there are not excitations and do not involve the physical 
magnetic field^'*.) This is further heightened by their being 
a rare instance of high-dimensional fractionalisation, a phe- 
nomenon of great interest in fields as diverse as correlated 
electrons and topological quantum computing^^. We hope our 
analysis will encourage our experimentalist colleagues to try 
their hand at directly detecting these monopole. There are 
surely many avenues to explore in search of useful signatures: 
scattering, transport and noise measurements, flux detection, 
to mention just a few. 
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APPENDIX A: THE DUMBELL PICTURE 

This material presents a detailed derivation of the dumb- 
ell Hamiltonian used extensively in our paper. We start from 
the generally accepted Hamiltonian which contains a sum of 
nearest-neighbour exchange and long range dipolar interac- 
tions. 



H = 



(y) L 



e^-Cj 3(e, •ry )(e^-T,^) 



SiSj . 
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The magnetic moment of a spin is denoted by fi, which 
equals approximately 10 Bohr magnetons (/i — lO/is) for 
the spin ice compounds discussed here (namely, Dy2Ti207 
and Ho2Ti207). The distance between spins is 7-y , and 
a ~ 3.54 A is the pyrochlore nearest-neighbour distance. 
D = noii^ / {Aira'^) = 1.41K is the coupling constant of the 
dipolar interaction. 

A dipole can be thought of as a pair of equal and opposite 
charges of strength ±q, separated by a distance a, such that 
/i = qd. The dipolar part of the Hamiltonian is reproduced 
exactly in the limit a ^ 0. Here we choose a to equal the 
diamond lattice constant aj, = \/3/2 a, therefore fixing q = 
ji/ad- The two ways of assigning charges reproduce the two 
orientations of the original dipole, as illustrated in Fig.|5] 

Given a configuration of N dipoles, let us label {g^, i = 
1, . . . , 2A^} the 2A^ charges in the corresponding dumbell 
configuration. The magnetic Coulomb interaction between the 
charges is given by 

ru 4 

47r Tij '-J ' 

voqiqj rij = 0, 



V(r, 



(A2) 



which reproduces the interaction between two dipoles exactly 
in the limit of large separation between them (modulo an in- 
nocuous '"self-energy" term discussed below). We can also 
tune the onsite contribution vq so as to match the interaction 
energy ±Jeff = ±(J + 5D)/3 between two neighbouring 
dipoles on the lattice (see for example Fig.|6]l. For two neigh- 
bouring spins both pointing into the shared tetrahedron (right 
panel of Fig. |6]l one has 



Jeff = w(0)-2w(ri2)-2w(r23)+v(ri3) 

= v{0) - v{ri2) - v{r23) + 2E^e\f + ^(ria), 



(A3) 
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FIG. 4: Phase diagram of spin ice in a [111] field. The location of the monopole liquid-gas transition from numerics (blue line) com- 
pared to experiment (black line). The solid line is a first-order transition terminated by a critical endpoint (fat circle). The dashed lines are 
crossovers. Inset: magnetisation curves showing the onset of first-order behaviour as the temperature is lowered. Our simulations cover 
the range 0.335 K < T < 0.8 K for 1024 spins. At the lowest temperatures, the parallel tempering code we use in our simulations of the 
Ewald-summed dipolar interaction no longer completely suppresses the hysteresis, and we have extended the first-order transition line using 
Clausius-Clapeyron. 



where v{r) — |V(r)|, and £^seif — —v{ri2) — —v{r23) is the 
self-energy of a dipole in the dumbell picture. For two neigh- 
bouring spins pointing into and out of the shared tetrahedron 
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FIG. 6: Schematic examples of two neighbouring dipoles on the 
pyrochlore lattice. 




(as in the left panel of Fig.|6]l, we obtain 

-Jeff = -v{0) + 2v{ri2) + 2v(r23) - v{ri3) 

= -v{0) + v{n2) + v{r23) + 2E,,ii - v{ri3). 

(A4) 



FIG. 5: Mapping from dipoles to dumbells. Left: two neighboring 
tetrahedra obeying the ice rule, with two spins pointing in and two 
out, giving zero net charge on each site. Right: The corresponding 
dumbell picture obtained by replacing each spin by a pair of opposite 
magnetic charges placed on the adjacent sites of the diamond lattice. 



Taking the difference between the two equations above, and 
using the fact that ri2 = ^23 and therefore v(ri2) = v{r23), 
one gets 



viO) = Jeff + 2u(ri2) -i;(ri3) 



(A5) 
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where 
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and in the end 



v{0) 



Mo 




2 

) 


■ 2 


1 " 


47r 








^ 2a 



(A6) 



Using the fact that D = /io/x^/47ra'^ and that a = ady/2/3, 
the expression above reduces to 



D. 



(A7) 



Now, the energy of a charge configuration can be straight- 
forwardly resummed and written in terms of the total charges 



Qa on each site a of the diamond lattice, i.e., Qa = Qii 
9i2 + 9i3 + for the quartet of charges with r^^ = r^^ 



The energy can thus be computed from a 



pairwise Coulomb interaction of the magnetic charges Qa'. 



47r Vc^is 



(A8) 



This gives the same total energy as Eq. IA2I up to an unim- 
portant overall constant (1/2) — ^if^/o-d)'^- Thus, it is 
also equivalent to the dipolar energy Eq. lAll up to corrections 
which are small everywhere, and vanish with distance at least 
as fast as l/r^ for each dipole pair. In particular, one can show 
that vq is large enough, for the compounds we consider, that 
it enforces Qa = everywhere in the ground state. These are 
just the ice rules, and the degeneracy of all the states follows 
from the fact that all terms with rap > vanish for Qa = 0. 
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